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■ - - Abstract 

ON , 

We investigate the problem of disseminating broadcast messages in wireless networks with time-varying links 
from a percolation-based perspective. Using a model of wireless networks based on random geometric graphs with 
dynamic on-off links, we show that the delay for disseminating broadcast information exhibits two behavioral 



(N 



■ regimes, corresponding to the phase transition of the underlying network connectivity. When the dynamic network 



is in the subcritical phase, ignoring propagation delays, the delay scales linearly with the Euclidean distance 
between the sender and the receiver. When the dynamic network is in the supercritical phase, the delay scales sub- 
linearly with the distance. Finally, we show that in the presence of a non-negligible propagation delay, the delay 
for information dissemination scales linearly with the Euclidean distance in both the subcritical and supercritical 
regimes, with the rates for the linear scaling being different in the two regimes. 



(N 



I. Introduction 

Large-scale wireless networks for the gathering, processing, and dissemination of information have 
become an important part of modern life. To ensure that important broadcast messages can be received 
^ ] by each node in a wireless network, the network needs to maintain full connectivity [1]. Here, the system 
• ; ensures that each pair of network nodes are connected by a path of consecutive links. In large-scale 

■ wireless networks exposed to severe natural hazards, enemy attacks, and resource depletion, however, 

O ■ 

. the full connectivity criterion may be overly restrictive or impossible to achieve. In these challenging 

^ ■ 

^ ! environments, the system designer may reasonably aim for a slightly weaker notion of connectivity, one 
c3 which ensures that a high fraction of the network nodes can successfully receive broadcast messages. This 
latter viewpoint can be explored using the mathematical theory of percolation [2]-[5]. 

In this paper, we investigate the problem of information dissemination in wireless networks from a 
percolation-based perspective. Using a model of wireless networks based on random geometric graphs 
with dynamic on-off links, we show that the delay for disseminating broadcast information exhibits a 
phase transition as a function of the underlying node density. Assuming zero propagation delay, we show 
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that in the subcritical regime, the delay scales linearly with the distance between the sender and receiver. 
In the supercritical regime, the delay scales sub-linearly with the distance. 

In recent years, percolation theory, especially continuum percolation theory [4], [5], has become a useful 
tool for the analysis of large-scale wireless networks [6]— [15]. A major focus of continuum percolation 
theory is the random geometric graph in which nodes are distributed according to a Poisson point process 
with constant density A, and two nodes share a link if they are within distance 1 of each other. A 
fundamental result of continuum percolation concerns a phase transition effect whereby the macroscopic 
behavior of the random geometric graph is very different for densities below and above the critical density 
A c . For A < A c (subcritical), the connected component containing the origin contains a finite number of 
points almost surely. For A > A c (supercritical), the connected component containing the origin contains 
an infinite number of points with a positive probability [3]-[5]. 

Wireless networks are subject to multi-user interference, fading, and noise. Thus, even when two nodes 
are within each other's transmission range, a viable communication link may not exist [7]. Furthermore, 
due to fading, the link quality can vary dynamically in time, inducing a frequently changing network 
topology. To capture these effects, we model a wireless network by a random geometric graph in which 
each link's functionality (activity) varies dynamically in time according to a Markov on-off process. Using 
this model, we investigate the problem of disseminating broadcast messages in wireless networks. Due to 
the dynamic on-off behavior of links, a delay is incurred in transmitting a broadcast message from the 
sender to the receiver even when propagation delay is ignored. The main question we address is how this 
delay scales with the distance between the sender and the receiver. 

As a first step, we show that the connectivity of the network with dynamic links exhibits a phase 
transition as a function of the underlying node density. We characterize the critical density for this phase 
transition in terms of the link state process. Next, we show that the delay for disseminating broadcast 
information exhibits two behavioral regimes, corresponding to the phase transition of the underlying 
network connectivity. When the dynamic network is in the subcritical phase, ignoring propagation delays, 
the delay scales linearly with the Euclidean distance between the sender and the receiver. This follows 
from the fact that in this regime, connectivity decays exponentially with distance, and on average, any 
information dissemination process is blocked by inactive links after the message travels a finite distance 
(and is resumed after the next link turns back on). When the dynamic network is in the supercritical 
phase, the delay scales sub-linearly with the distance between the sender and the receiver. In this case, the 
delay is determined largely by the amount of time it takes for the message to reach the infinite connected 
component of the dynamic network. Finally, we characterize the delay for information dissemination when 
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propagation delays are taken into account. Here, the problem becomes more subtle. We show that, with 
the presence of a non-negligible propagation delay, the delay for information dissemination scales linearly 
with the Euclidean distance between the sender and the receiver in both the subcritical and supercritical 
regimes, with the rates for the linear scaling being different in the two regimes. 

In order to study the behavior of information dissemination delay in wireless networks with dynamic 
links, we model the problem as a first passage percolation process [16], [17]. Similar first passage 
percolation problems have been studied within the context of lattices [3], [16]. Related continuum models 
are considered in [8], [13], [17]. In [17], Deijfen studies a continuum growth model for a spreading 
infection with Poisson point processes, and shows that the shape of the infected cluster scales linearly 
with time in all directions. In [8], Dousse et al. study how the latency of information dissemination 
scales within an independent site percolation model in wireless sensor networks. There, each sensor 
independently switches between the on and off states at random from time to time. The authors show 
that the latency scales linearly with the distance between the sender and the receiver when the dynamic 
sensor network is in the subcritical phase. In [13], the authors obtain similar results for degree-dependent 
site percolation model in wireless sensor networks. Unlike the problems studied in [8], [13], however, the 
problem addressed in this paper requires a bond percolation model, which demands different modelling 
and analysis techniques. Furthermore, in contrast to [8], [17], we also study the delay scaling for networks 
in the supercritical phase. Finally, we present new results regarding networks with propagation delay. 

The remainder of this paper is organized as follows. In Section II, we outline some preliminary results 
for random geometric graphs and continuum percolation. In Section III, we present a simple model for 
wireless networks with static unreliable links. In Section IV, we introduce a more sophisticated model for 
wireless networks with dynamic unreliable links, and present our main results regarding percolation-based 
connectivity and information dissemination within this model. In Section V, we present simulation results, 
and finally, in Section VI, we conclude the paper. 

II. Random Geometric Graphs and Continuum Percolation 
A. Random Geometric Graphs 

We use random geometric graphs to model wireless networks. That is, we assume that the network 
nodes are randomly placed over some area or volume, and a communication link exists between two 
(randomly placed) nodes if the distance between them is sufficiently small, so that the received power 
is large enough for successful decoding. A mathematical model for this is as follows. Let || ■ || be the 
Euclidean norm, and /(•) be some probability density function (p.d.f.) on IR d . Let X 1 ,X 2 ,...,X„ be 
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independent and identically distributed (i.i.d.) o?-dimensional random variables with common density /(•), 
where Xj denotes the random location of node i in M. d . The ensemble of graphs with undirected links 
connecting all those pairs {x*, x,-} with ||xj — Xj|| < r,r > 0, is called a random geometric graph [5], 
denoted by G(X n) r). The parameter r is called the characteristic radius. 

In the following, we consider random geometric graphs G(X n , r) in M 2 , with X 1; X 2 , X n distributed 
i.i.d. according to a uniform distribution in a square area A = [0, a/^] 2 . Let A = \A\ be the area of A. 
There exists a link between two nodes i and j if and only if i lies within a circle of radius r around Xj. 
As n and A both become large with the ratio ^ = A kept constant, G(X n ,r) converges in distribution 
to an (infinite) random geometric graph G(H\,r) induced by a homogeneous Poisson point process with 
density A > 0. Due to the scaling property of random geometric graphs [4], [5], we focus on G(H\, 1) 
in the following. 

B. Critical Density for Continuum Percolation 

To intuitively understand percolation processes in large-scale wireless networks, consider the following 
example. Suppose a set of nodes are uniformly and independently distributed at random over an area. All 
nodes have the same transmission radius, and two nodes within a transmission radius of each other are 
assumed to communicate directly. At first, the nodes are distributed according to a very small density. This 
results in isolation and no communication among nodes. As the density increases, some clusters in which 
nodes can communicate with one another directly or indirectly (via multi-hop relay) emerge, though the 
sizes of these clusters are still small compared to the whole network. As the density continues to increase, 
at some critical point a huge cluster containing a large portion of the network forms. This phenomenon 
of a sudden and drastic change in the global structure is called a phase transition. The density at which 
phase transition takes place is called the critical density [3]-[5]. 

More formally, let H\ : o = H\ U {0}, i.e., the union of the origin and the infinite homogeneous Poisson 
point process with density A. Note that in a random geometric graph induced by a homogeneous Poisson 
point process, the choice of the origin can be arbitrary. We have the following definition [4]. 

Definition 1: For G(H\,o : 1), let W be the connected component of G(H\ : o, 1) containing 0. Define 



SUBMITTED TO IEEE TRANSACTIONS ON INFORMATION THEORY 



5 



the following critical densities: 

A # 4 inf{A : Pr(|W„| = oo) > 0}, (1) 

\ N 4 inf{A : E[\W \] = oo}, (2) 

A c = inf{A : Pr(d(W ) = oo) > 0}, (3) 

X D 4 inf{A : £[d(W )] = oo}, (4) 

where |W | is the cardinality — the number of nodes — of Wq, and d(W ) = sup{||x — y|| : x, y e Wo}. 

As shown in Theorem 3.4 and Theorem 3.5 in [4], these four critical densities are identical. According 
to the theory of continuum percolation [4], < A c < oo. Furthermore, when A > A c , there exists a unique 
infinite component in G(H\ : o, 1) with probability 1, and when A < A c , there is no infinite component in 
G(Hx, , 1) with probability 1 [4]. 

III. Wireless Networks with Static Unreliable Links 

Random geometric graphs are good simplified models for wireless networks. However, due to noise, 
fading, and interference, wireless communication links between two nodes are usually unreliable. We first 
use the bond percolation model on random geometric graphs to study percolation-based connectivity of 
large-scale wireless networks with static unreliable links. Given a random geometric graph G(H\, 1), let 
each link of G(H\, 1) be active (independent of all other links) with probability p e (d) which may depend 
on d, where d = ||xj — Xj|| < 1 is the length of the link The resulting graph consisting of all active 
links and their end nodes is denoted by G(H\, l,p e (-)). This model is a specific example of the random 
connection model in continuum percolation theory [4]. In this simple model, all links in the network are 
either active (on) or inactive (off) for all time. Later in this paper, we will study a more sophisticated 
model where links dynamically switch between active and inactive states from time to time. 

Definition 2: For G(H\ : o, l,p e (-)), let Wq be the connected component of G(H\,o, l,p e (-)) containing 



0. We define four critical densities: 

A # (p e (0) = inf{A : PrflWft = oo) > 0}, (5) 

A JV (p e (.)) = inf{A : E[\W^\] = oo}, (6) 

A c (p e (.)) = inf{A : Pv(d(W^) = oo) > 0}, (7) 

X D ( Pe (-)) ± inf{A : E[d{W' )\ = oo}, (8) 



where \Wq\ is the cardinality — the number of nodes — of Wq, and d(W$) = sup{||x — y|| : x, y e Wq}. 
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As in traditional continuum percolation, the following proposition asserts that the above four critical 
densities are identical. 

Proposition 1: For G(Hx,o, l,Pe( - ))' we nave 

A # (p e (0) = Aiv(Pe(-)) = A c (Pe(0) = ^(Pe(0)- (9) 

Proof: The identity A # (p e (-)) = A^(p e (-)) is given by Theorem 6.2 in [4]. 

We now show A#(p e (-)) = A c (p e (-)). The proof method is similar to the one used for Theorem 3.4 
in [4]. Suppose A > A # (p e (-)). Then for some 5 > 0, Pr(|W^| = oo) = 5 > 0. For every h > 0, the 
box = [— /i, /i] 2 contains at most a finite number of nodes of G(H\ : o,l,p e (-)) with probability 1. 

Thus, Pr(|W£ n 5(/i) c | = cx)) = 5 > 0. However, {\W^ n 5(/i) c | = oo} implies {|W£ n 5(/i) c | > 0}, 
so that d(Wo) > h. Hence we have Pr(d(Wo) > /i) = 5 > 0. Since this holds for all h > 0, we have 
A > A c (p e (-))- Therefore, A # (p e (0) ^ A c (p e (-)). 

To show A # (p e (-)) < A c (p e (-)), note that cL(Wq) < \Wq\ — I, where equality is obtained when Wq is a 
chain and the distance between any two adjacent nodes equals 1. Thus, {\Wq\ < oo} implies {cL(Wq) < 
oo}. This proves A # (p e (0) = A c (p e (-))- 

Finally, we show A D (p e (-)) = Ajv(p e (-)). Since < |W^| - 1, {E[d(W^)] = oo} implies 

{-EOWqI] = oo}. Thus we have X D (p e (-)) > X N (p e (-)). On the other hand, if A > X N (p e (-)), then 
A > A c (p e (-)), i.e., Pr(d(WQ) = oo) > 0. As a consequence, E[d(WQ)} = oo, which implies Ajv(p e ( - )) > 
A D (p e (-))- Therefore, A D (p e (-)) = Aiv(p e (0)- D 

Since the four critical densities are identical, in the remainder of this paper, we state our results with 
respect to A c (p e (-)). 

It is known that when A > A c (p e (-)), G(H\, l,p e (-)) is percolated, i.e. with probability 1, there exists 
a unique infinite component in G(H\, 1) consisting of active links and their end nodes, and when A < 
A c (p e (-))' G(H\, l,p e (-)) is not percolated, i.e., with probability 1, there is no infinite component in 
G(H\, 1) consisting of active links and their end nodes [4]. 

The following monotonic property for A c (p e (-)) can be easily proved by coupling methods. 

Proposition 2: Let A c (p e (-)) and X c (p' e (-)) be the critical densities for G{H\ 1 l,p e (-)) and G(H\, l,p' e (-)), 
respectively. Then, if p' e (x) < p e (x),yx e (0, 1], we have A c (p e (-)) < A c (p' e (-)). 

The following proposition asserts that when the random connection model is in the subcritical phase, 
the probability that the origin and a given node are connected decays exponentially with the distance 
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between them. This is analogous to similar results in traditional continuum percolation (Theorem 2.4 
in [4]) and discrete percolation (Theorem 5.4 in [3]). 

Proposition 3: Given G(H x ,oA,Pe(-)) with A < A c (p e (-)), let B(h) = [-h,hf, h e R + . Then there 
exist constants ci,c 2 > 0, such that Pr(0 B(h) c ) < c\e~ C2h , where {0 B(h) c } denotes the event 
that the origin and some node in B(h) c are connected, i.e., the origin and some node outside B(h) are 
in the same component. 

The proof for this proposition is similar to the one for Theorem 2.4 in [4]. For completeness, we give 
the proof in Appendix A. 



A. Percolation-based Connectivity 

For the random connection model, we assumed that the structure of the graph does not change with 
time. Once a link is active, it remains active forever. In wireless networks, however, the link quality 
usually varies with time due to shadowing and multi-path fading. In order to study percolation-based 
connectivity of wireless networks with time-varying links, we investigate a more sophisticated model. 
Formally, given a wireless network modelled by G(H\, 1), we associate a stationary on-off state process 
{Wij(dij,t);t > 0} with each link where dij is the length of the link, such that Wij(dij,t) = if 

link (i, j) is inactive at time t, and Wij(dij,t) = 1 if link is active at time t. A similar problem for 
discrete lattice has been studied in [18]. Our model can be viewed as one of dynamic bond percolation 
in random geometric graphs. 

For such dynamic networks, we will show that there exists a phase transition, and the critical density 
for this model is the same as the one for static networks with the corresponding parameters. To simplify 
matters, assume that {Wij(dij,t)} is probabilistically identical for all links with the same length. Use 
{W(d, t)} to denote the process for a link with length d when no ambiguity arises. Assume that {W(d, t)} 
is a Markov on-off process with i.i.d. inactive periods Y k (d), k > 1, and i.i.d. active periods Z k (d), k > 1, 
where E[Y k (d) + Z k (d)] < oo, Vi(Z k (d) > 0) = 1 and Pr (Y k (d) > 0) = 1 for < d < 1. That is, 
both the active and inactive periods are always nonzero. Further assume that mf 0<d <i{E[Y k (d)}} > and 
^Po<d<i{E[ Y k(d)}} < oo. 

Under the above assumptions, the stationary distribution of {W(d,t)} is given by [19] 



IV. Wireless Networks with Dynamic Unreliable Links 





(10) 



Vo (d) ±Pr(W(d,t) =0) 



E[Z k (d)] + E[Y k (d)Y 
E[Y k (d)} 



(ID 



E[Z k (d)}+E[Y k (d)Y 
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where 771(d) is the active ratio for a link with length d. 

Let the graph at time t be G(TC\, 1, W(d, t)). That is, G(TCx,l,W(d,t)) consists of all active links 
at time t, along with their associated end nodes. The following theorem establishes a phase transition 
phenomenon with respect to connectivity in a wireless network with dynamic unreliable links modelled 
by G(TC\, 1, W(d,t)). It also asserts that the critical density is the same as the one for the static network 
G(H\, 1, 771(d)), i.e, the network in which each link is active with probability 771(d). 

Theorem 4: Let X c (r]i(d)) be the critical density for the static model G(H\, 1, 771(d)). Then G(H\, 1, W(d, t)) 
is percolated for all t > if A > A c (771(d)), and not percolated at any t > if A < A c (77i(d)). 

Proof: Since A > A c (77i(d)) and < 771(d) < 1, Vd G (0, 1], by the monotonic property of X c (p e (-)) 
(Proposition [2]), we can construct a new model G(H\,l,W'(d,t)) and choose e > such that A > 
A c (7?i(d)) > A c (77i(d)) and < n[(d) < 1, Vd G (0, 1], where ^(d) = (1 - 6)771 (d), for d e (0, 1]. As 
active periods are always nonzero, we can choose 5 > such that for any link 

Pr(W ij (S) = l\W ij (d,0) = l)>l-e, 

where Wij(8) = min tg [ 0i< 5] Wij(d,t). Then, 

Pr(^(5) = l)> (1-6)7/1 (d) = V [(d). 

Since A > A c (?7i(d)), for any t E [0,5], G(H\,l,W(d,t)) is percolated. Repeat this argument for all 
intervals [kS, (k + 1)5} with integer k. Let E k be the event that G(H\, 1, W(d, t)) is percolated for all 
t E [k5, (k + 1)5]. Then, we have 

Prffl^J =l-Prf|J^ C ) >l-E Pr ^) = 1 - 

Similarly, when A < A c (?7i(d)), we can construct another model G(H\, 1, W"(d,t)) and choose e > 
such that A < X c (rj'{(d)) < K(vi(d)) and < r)'{(d) < l,Vd e (0,1], where 7]'{{d) = e(l - 771(d)) + 
771(d), Vd G (0, 1]. Since inactive periods are always nonzero, we can choose 5 > such that for any link 

Pr(W ii (5) / = 0|Wy(d,0) = 0) > 1 - e, 

where Wy(5)' = max 4g [ 0j(5 ] Wij(d, t). Then, 

Pr(W^)' = 0) < 1 - (1 - 77i(d))(l - e) = 77i'(d). 

Since A < \ c (r)1(d)), for any t G [0, 5], G(H\, l,W(d,t)) is not percolated. Repeat this argument for all 
intervals [fc<J, (A; + 1)5] with integer k, and then proceed in the same way as before, i.e., using countable 
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When the process {W(d, t)} is independent of link length d, we use {^(t)} to denote the process, 
and 771 and r] to denote its stationary distribution. 

B. Information Dissemination in Wireless Networks with Dynamic Unreliable Links 

We have shown that there exists a critical density X c (rji(d)) such that when A > \ c (Vi(d)), G(H\, 1, W(d, t)) 
is percolated for all time. If G(H\, l,W(d,t)) is percolated, when one node inside the infinite component 
of G(TC\, 1, W(d, t)) broadcasts a message to the whole network, then assuming that there is no propaga- 
tion delay, all nodes in the infinite component of G(H\, 1, W(d ) t)) receive this message instantaneously. 
The nodes in the infinite component of G(H\, 1) but not in the infinite component of G(H\, 1, W(d, t)) 
cannot receive this message instantaneously. Nevertheless, as links switch between the active and inactive 
states from time to time, those nodes can still receive the message via multi-hop relaying at some later 
time. This remains true even if A < \ c (r]i(d)) and G(H\, 1, W(d, t)) is never percolated. In this case, 
when one node inside the infinite component of G(H\, 1, W(d, t)) broadcasts a message, due to poor 
connectivity, only a small number of nodes can receive this message instantaneously. However, as long as 
two nodes u and v are in the infinite component of G(H\, 1), the message can eventually be transmitted 
from u to v over multi-hop relays. The main question we address here is the nature of this information 
dissemination delay. 

This problem is similar to the first passage percolation problem in lattices [3], [16]. Related continuum 
models are considered in [8], [13], [17]. In [17], the author study continuum growth model for a spreading 
infection. In [8] and [13], the authors consider wireless sensor networks where each sensor has independent 
or degree-dependent dynamic behavior, which can be modelled by an independent or a degree-dependent 
dynamic site percolation on random geometric graphs, respectively. The main tool is the Subadditive 
Ergodic Theorem [20]. We will use this technique to analyze our problem. 

In the following, we will show that in a large-scale wireless network with dynamic unreliable links, 
the message delay scales linearly with the Euclidean distance between the sender and the receiver if the 
resulting network is in the subcritical phase, and the delay scales sub-linearly with the distance if the 
resulting network is in the supercritical phase. 

To begin, we define the delay on a link as the amount of time for node i to deliver a packet 
to node j over link (i, j). In particular, ignoring propagation delay, if is active when % initiates a 
transmission, then the delay is zero. If is inactive, the delay is the time from the instant when % 
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initiates transmission until the instant when becomes active. Mathematically, let delay Zy(c?y) be a 

random variable associated with link having length dij, such that 

Pr(2y(dy) = 0) = 771 (dy), 9 
Pr(T y (dy)>t) = ^KOP^W, UZj 

where P dij (t) = Pr(Wy(dy, f) = 0,Vtf G [0, t)| (t%, 0) = 0), and (771(d), 770(d)) is the stationary 

distribution of {W(d,t)} given by CUS) and (fTTT) . 

Let d(u, v) = | |X„ — X„| | and 

T(u, v) = T(X U , X„) 4 inf J V T ii (^) I , (13) 

where l(u,v) is a path of adjacent links from node u to node v, and C(u,v) is the set of all such paths. 
Hence, T(u, v) is the message delay on the path from u to v with the smallest delay[] 

Theorem 5: Given G(H\,l,W(d,t)) with A > A c , there exists a constant 7 satisfying 7 < 00 and 
7 > with probability 1, such that for any u, v G C(G(H\, 1)), where C(G(H\, 1)) denotes the infinite 
component of G(H\, 1), 

(i) if G(Ti\, 1, W(d, £)) is in the subcritical phase, i.e., A < \Jj}\{d)), then for any e > 0, 5 > 0, there 
exists d < 00 such that for any w, u with d(u, v) > do, 

T(u,v) 



Pr 



7 



< e > 1 - 8; (14) 



d(tt, v) 

(ii) if G(H\, 1, W(d, i)) is in the supercritical phase, i.e., A > X c (rji(d)), then for any e > 0, 5 > 0, there 
exists d < 00 such that for any u, v with d(u,v) > d , 

Pr(^4<e)>l-5. (15) 



d(u, v 

Before proceeding, we introduce some new notation. Let 

X t 4 argmin {| | Xj - (i, 0)| |}, (16) 

X je C(G(W A ,l)) 

T z , m 4 T(X,X m ), for ||X-X m || < oo,0 < / < m. (17) 
The proof for Theorem [5]-(i) is based on the following lemma: 
Lemma 6: Let 

A lim (18) 

m— *oo 777, 



Then, 7 = inf m >i E[T °; m] , and lim^^ ^ = 7 with probability 1. 



'Note that the path with the smallest delay may be different from the shortest path (in terms of number of links) from node u to node v. 
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To show Lemma [61 we use the following Subadditive Ergodic Theorem by Liggett [20]. 

Theorem 7 (Liggett [20]): Let {5j, m } be a collection of random variables indexed by integers < / < 
m. Suppose {5/, m } has the following properties: 

(i) So, m < S 0; i + S lttn , < I < m; 

(ii) {Sf m -i)k,mk> m > 1} is a stationary process for each k; 

(iii) {Sij + k, k > 0} = {Si+ij+k+i, k > 0} in distribution for each I; 

(iv) E'fjS'cml] < oo for each m. 
Then 

(a) a 4 hm m ^ 00 ^1 = inf m >! S 4 lim^ %p exists with probability 1 and E[S] = a. 
Furthermore, if 

(v) the stationary process in (ii) is ergodic, 
then 

(b) S = a with probability 1. 

To show Lemma [6l we need to verify that the sequence {T/ m , I < m} satisfies conditions (i)-(v) of 
Theorem |7J It is easy to see that (i) is satisfied, since T m is the delay of the path with the smallest 
delay from X to X m and T 0j / + T^ m is the delay on a particular path from X to X/ (it has the smallest 
delay from X to X/, and from X; to X m ). Furthermore, because all nodes are distributed according to 
a homogeneous Poisson point process, the geometric structure is stationary and hence (ii) and (iii) are 
guaranteed. We need only to show conditions (iv) and (v) also hold for {TJ jm , I < m}. To accomplish 
this, we first show property (iv) holds for {T i m ,/ < m}. 

Lemma 8: Let r = | |X — (0, 0)| |, then r < oo with probability 1. 

Proof: We consider a mapping between G(H\, 1) and a square lattice C = d ■ Z 2 , where d is the edge 
length. The vertices of C are located at (d x i, d x j) where 6 Z 2 . For each horizontal edge a, let 
the two end vertices be (d x a x ,d x a y ) and (d x a x + d, d x a y ). 

For edge a in C, define event A a (d) as the set of outcomes for which the following condition holds: 
the rectangle R a = [a x d — |, a x d + ™] x [a y d — |, a y d + |] is crossec 2 ] from left to right by a connected 
component in G(H\, 1). If A a (d) occurs, we say that rectangle R a is a good rectangle, and edge a is a 

2 Here, a rectangle R = [xi, X2] x [2/1,2/2] being crossed from left to right by a connected component in G(7i\, 1) means that there exists a sequence of 
nodes vi, t>2, v m G G(H\,1) contained in R, with \\x Vi — x lli+1 || < 1, i = l,...,m — 1, and < x(vi) — x\ < 1,0 < X2 — x(v m ) < 1, where 
x(vi) and x(v m ) are the i-coordinates of nodes v\ and v m , respectively. A rectangle being crossed from top to bottom is defined analogously. 
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(a) Good Rectangle (b) Open Rectangle 

Fig. 1. Examples of good and open rectangles (edges) 

good edge. Let 

p g (d) ± Pr(A a (d)). 

Define A a (d) similarly for all vertical edges by rotating the rectangle by 90°. An example of a good 
rectangle and a good edge is illustrated in Figure [T]- (a). 

Further define event B a (d) for edge a in £ as the set of outcomes for which both of the following 
hold: (i) A a (d) occurs; (ii) the left square S~ = [a x d — |, a x d + |] x [a y d — |, a y d + |] and the right 
square S+ = [a x d + ^, a x d + ^] x [a y d — |, a y d + |] are both crossed from top to bottom by connected 
components in Gi(H\, 1). 

If B a (d) occurs, we say that rectangle R a is an open rectangle, and edge a is an open edge. Let 

Po (d) ^Pr (B a {d)). 

Define B a (d) similarly for all vertical edges by rotating the rectangle by 90°. Examples of an open 
rectangle and an open edge are illustrated in Figure [U-(b). 

Suppose edges b and c are vertically adjacent to edge a, then it is clear that if events A a (d), A b (d) and 
A c (d) all occur, then event B a (d) occurs. Moreover, since events A a (d), At(d) and A c (d) are increasing 
eventfl by the FKG inequality [3]-[5], 

Po{d) = Pv(B a (d)) 

> Pr(A a (d) n A b (d) n A c (d)) 

> Px{A a {d))Px(A b {d))Pi(A c (d)) 

= (pM) 3 - 

According to Corollary 4.1 in [4], the probability p g (d) converges to 1 as d — > oo when G(H\, 1) is in 
the supercritical phase. In this case, (p g (d)) 3 converges to 1 as d — >• oo as well. Hence, p (d) converges 
to 1 as d — >• oo when G(H\, 1) is in the supercritical phase. 

3 An event A is called increasing if Ia(G) < Ia(G') whenever graph G is a subgraph of G' , where Ia is the indicator function of A. An event A is 
called decreasing if A c is increasing. For details, please see [3]-[5]. 
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Note that in our model, the events {B a (d)} are not independent in general. However, if two edges 
a and b are not adjacent, i.e., they do not share any common end vertices, then B a (d) and Bb(d) are 
independent. Furthermore, when edges a and b are adjacent, B a (d) and B b (d) are increasing events and 
thus positively correlated Q. Consequently, our model is a 1 -dependent bond percolation model. It is known 
that there exists P^°dep < ^ sucn mat an y 1 -dependent model with p > p\™ep * s P erc °l ate d, where p is 
the probability of an edge being open [21]. 

Now define 

d 4 inf jd' > 1 : p {d') > max |^# d n e d p || , (19) 

and choose the edge length of C to be d > d . Then there is an infinite cluster consisting of open edges 
and their end vertices in C. Denote this infinite cluster by C(C). 

From Figure [2] it is easy to see that all the nodes along the crossings in R a and all the nodes along the 
crossings in Rb for any a, b G C(C) are connected. Since the infinite component of G(H\, 1) is unique, 
all the nodes along the crossings in R a for each a G C(C) must belong to C(G(H\, 1)). 

By definition, no node of G(TC\, 1) strictly inside A(0,r ) belongs to C(G(TC\, 1)). This implies that 
no edge of C strictly inside .4.(0, r ) belongs to C(C). To see this, suppose edge a^j of C is strictly inside 
4.(0, r ) and belongs to C(C). The nodes along the crossings in R a .. belong to C(G(TC\, 1)). As shown 
in Figure (3]-(a), when d > 1 and r ^> 1, no matter what direction the edge a^j has, there are some nodes 
along the crossings in R a . . (therefore belonging to C(G(H\, 1))) which are strictly inside 4.(0, r ). These 
nodes then have strictly smaller distance to than node X . This contradiction ensures that no edge of 
C strictly inside 4.(0, r ) belongs to C(C). 

"Positive correlation means Pi(B a (d)\B b (d)) > Pr(B a (d)). 
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Consider the dual lattice £' of C. The construction of £ is as follows: let each vertex of CJ be located 
at the center of a square of C. Let each edge of £' be open if and only if it crosses an open edge of C, 
and closed otherwise. It is clear that each edge in C is open also with probability p (d). Let 

q = l-p {d) < i. 

Choose 2m edges in C . Since the states (open or closed) of any set of non-adjacent edges are 
independent, we can choose m edges among the 2m edges such that their states are independent. As 
a result, 

Pr(all the 2m edges are closed) < q m . 

Now a key observation is that if no edge of C strictly inside A(0,r ) belongs to C(C), for which 
the event is denoted by Ec, then there must exist a closed circuit in C (a circuit consisting of closed 
edges) containing all edges of C strictly inside A(0,r ), for which the event is denoted by Ec and vice 
versa [3]. This is demonstrated in Figure |3]-(b). Hence 

Pr(E c ) = 1 ^ Pt{E £ ,) = 1. 

Any closed circuit in £' containing all edges of C strictly inside ^4(0, r ) has length greater than or 
equal to 21, where I = 2 [^J . Thus we have 

oo oo 

Pi(E £t ) = ^Pr(3C c (2m)) < ^ 7 (2m)<f\ 

m=l m=l 
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where O c (2m) is a closed circuit having length 2m in £ containing all edges of C strictly inside .4.(0, r ), 
and 7 (2m) is the number of such circuits. By Proposition [T51 in Appendix B, we have 7 (2m) = 



m — 



l)3 2m so that 



X>(2m)g m < £_(m-l)(9ff) 



m=l 



m=l 

4 

27 
4[i 



m=l 



m=l 



(/ - 2)9g] 



(9g)'. 



(20) 



27(1 -9g) 2 

Since g < |, we have Pt(Ec) — > as / = 2L^J — >■ 00. That is, as r goes to infinity, with probability 
1, there is some edge of C strictly inside 4.(0, r ) belonging to C(C). Hence, with probability 1, there is 
some node of G(H\, 1) strictly inside 4.(0, r ) belonging to C(G(H\, 1)). This contradiction implies that 
r is finite with probability 1. □ 



Let r m = ||X m — (m, 0)||, by Lemma [8] and stationarity, we have r m < 00 with probability 1, for any 
m. 

Lemma 9: Let L(X , X m ) be the shortest path (in terms of the number of links) from X to X m , and 
let |L(X , X m ) I denote the number of links on such a path. If | |X — X m | | < 00, then |L(X , X m ) | < 00, 
and E[T^ m ] < 00, where T L m denotes the delay on path L(X ,X m ). 

Proof: We use the same mapping as the one for the proof of Lemma [U For any given ^/| < 5 < 1, 
define 

d s = max{M{d' : p g (d') > 6},\\X - X m ||}. (21) 

Then, for any d > ds, we have p g (d) > 5. 

Now, consider a fractal structure as shown in Figure @] first a square S(ds) is constructed with edge 
length dg centered at X °+ Xm . Then, a second square S(3ds) is constructed with edge length 3^5 also 
centered at X °+ Xm . The construction proceeds in the same manner, i.e., at step j, a square S(3-'^ 1 ds) is 
constructed with edge length &~ l ds centered at x °^ Xm . Thus, we have the initial square and a sequence 
of square annuli that do not overlap. 

Denote the square annulus with inside edge length 3^ 1 ds (j > 2) and outside edge length 3^ds by 
D(3 j d s ). Let Aj be the event that the upper horizontal rectangle of D(3 j d s )— [f - f d s , f + %d s ] x 
[^2~c?5, y^s] is good, i.e., it is crossed by a connected component in G(H\, 1) from left to right. Since 
the length of the corresponding lattice edge of the upper horizontal rectangle of D(3^ds) is 2-3^~ 1 ds > ds, 
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Fig. 4. Square annuli 

we have Pr{A^} > 5. Similarly define Aj, B^~ and Bj to be the events that the lower, right and left 
rectangles are good, respectively. Then Pr{Aj} > 5, Pr{£> + } > 5 and Vr{Bj} > 5, Vj > 1. 

Let Ej be the event that there exists a circuit of connected nodes in G(TC\, 1) within Dffids). Once 
Aj,Aj,Bj~ and Bj all occur, Ej must also occur. Although A^j,Aj,Blj and Bj are not independent, 
they are increasing events. By the FKG inequality, we have 

Px(Ej) > Pr(A-~ n Aj n Bf n Bj) 

> Pr(A+) Pr (Aj) Pr(Bf) Pr (Bj) 

> 5 A . (22) 

When Ej occurs, X and X m are contained in S(3 :j ~ 1 ds) and there is a circuit of connected nodes in 
G(TC\, 1) contained in the square annulus Dftds). If the shortest path between X and X m , L(X , X m ), 
were to go outside S(3^dg), it would intersect the closed circuit contained by D^dj) and we could 
construct a shorter path from X to X m . This implies that L(X ,X m ) must be contained in S(?>id$). 

Suppose u, v and w are three consecutive nodes along L(X ,X m ). Then ||X U — X^H > 1, since 
otherwise v would not belong to the shortest path. Hence, if we draw circles with radius |, centered at 
X M and X w , respectively, then the two circles do not overlap. Consequently, if the length of L(X , X m ) 
is \L\ = |L(X ,X m )|, then we must be able to draw at least [^-] circles with radius | centered at 
alternating nodes along L(X ,X m ). All of these circles are contained in the square with edge length 
3 j ds + 1. Such a square contains at most \(3 j ds + l) 2 /[7r(|) 2 ]] non-overlapping circles with radius |. 
Therefore, \L\ < 2\4(3 j d s + 1) 2 /tt] < oo. 

Now if \L\ > 2\4(3 j d 5 + l) 2 /n], then \L\ > 2\4(3 l d 5 + 1) 2 /tt] for all i = 1, 2, j. By the above 
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argument, none of the events E u E 2 , ...Ej can occur. Thus 

) <f[rrm<(i-5y. 

' i=l 



Pr ML| > 2 
Let M = 2 ["±(3d* + l) 2 ] , then we have 

oo 

E[\L\] = ^Pr(|L|>fc) 



(&d s + If 



k=0 

M 



^Pr(|L| > k) + Pr d L l > k ) 



k=0 



k=M+l 



< 



j=i 

oo 



7T 



•(••V 'r/, • I)'- 



Pr \L\ > 



7T 



: (3H + i; 



/a \ 

< M + J^l -(3^ + l) 2 + lj (i-sy 

j=i ^ ' 

= M + J2 ( -(9 ■ + 6 • 3 J 'd 5 + !) + ! (1- 
j=i ^ ' 

= M + ^^(i~6y + ^^(i-5y+(^ + i\f^(i-5y. (23) 

71 .7 = 1 71 .7 = 1 ^ ' .7 = 1 



OO. 



When 5 > ^/§, we have (1 - t^ 4 )-? < 9^'. Thus, < 
Let A^) = sup 0<d < 1 {?7o(rf) J B[Y fe (c/)]} < oo, then 

E[T L m \\L\] = J2v$\d)E[Y®(d)} < \L\A wm , (24) 

i=l 

where ^(d) and E[Y^\d)] are the stationary probability of the inactive state, and the expected inactive 
period of the ?-th link with length d on L(X , X m ), respectively. Hence 

E[T L J = E[E[T L m \\L\]] < E[\L\]A wm < oo. (25) 

□ 

To show property (v), we show {T( jri ~i)j,mj ) m > 1} is strong mixing^ 

Lemma 10: The sequence {T , ( m _ 1 ) fc mfc , m > 1} is strong mixing, so that it is ergodic. 

Proof: From the proof of Lemma [81 we have Pr(-E'j) > <5 4 for all j = 1, 2, .... Summing over j yields 

oo oo 

^Pr(^)>^5 4 = oo. (26) 

0=1 3=1 



A measure preserving transformation H on (f2, P) is called strong mixing if for all measurable sets A and B, linim—joo \P(A f~l H~ m B) — 
P(A)P(B)\ = 0. A sequence {X n ,n > 0} is called strong mixing if the shift on sequence space is strong (weak) mixing. Every strong mixing system is 
ergodic [22]. 
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Fig. 5. As k — > oo, the paths inside A\ and A2 do not share any common nodes. Hence T( jn _x)j, r aj an d T( m +k-i)jJm+k)j are independent 
of each other as k — » 00. 

Since £j are independent events, by the Borel-Cantelli Lemma, with probability 1, there exists f < oo 
such that Eji occurs. 

We now construct squares A l and A 2 centered at X(m " 1) 2 J+Xmj and X(m+fc - 1)j 2 +X(m+fc)j with edge length 
3i'ds and 3 J "ds respectively, such that the path with the smallest delay from X( m _x)j to X m j, and the 
path with the smallest delay from X( TO+ fc_i)j to X( m+fc )j are contained in At and A 2 , respectively. Let E 
be the event that f < 00 and j" < 00. Then Pt(E) = 1. 

When finite j' and j" exist, due to stationarity, j' and j" are independent of k. Hence, as k — > 00, 
Ai and A 2 become non-overlapping so that the paths inside A\ and A 2 do not share any common nodes 
of G(H\, 1). Hence T( m -i)j, m j an d ^( m +fc-i)j,(m+fc)j are independent of each other as k — > 00. This is 
illustrated in Figure [51 

Therefore 

lim Pr({T (m _ 1)ijmi < t} n {T ( 
= lim Pr({T (m _ 1)ijmi < t} H {T (m+fc _ 1)ji(m+A . )i < Pr(£) 

fc^oo 

+ lim Pr({T (m _ 1)i>mi < t} n {T (m+fe _ 1)i)(m+fe)j < i'}|£ c ) Pi{E c ) 

< t'\E) 

= P^(T( m -i)j jm j < t) Pr(T( m _!)j im j < i), (27) 
This implies that sequence {T( m _i) fc)mfc , m > 1} is strong mixing, so that it is ergodic. □ 

Now, we present the proof for Lemma [6] 

Proof of Lemma [6J - Conditions (i)-(iii) of Theorem [7] have been verified. The validation of (iv) is 
provided by Lemma |9] Let L(X , X m ) be the shortest path from X to X m . Since L(X , X m ) is a 
particular path, we have T 0im < T^ m so that E[T 0>m ] < E[T^ m ], where T^ m denotes the delay on path 
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Fig. 6. Path segments of the paths from Xo to X m . 

L(X , X m ). By Lemma[9l we have E[TQ m ] < oo and therefore E[T 0tm ] < oo. Furthermore, due to Lemma 
[TOl {T( m _!) fe)mfc , m > 1} is ergodic, thus the results (a) and (b) of Theorem [7] hold. □ 

Remark: Using the proof for condition (iv) of Theorem [71 we can show that for any two nodes u and 
v in the infinite component of G(H\, 1) which are within finite Euclidean distance of each other, i.e., 
m,d6 C(G(H\, 1)) with d(u,v) < oo, E[T(u,v)] < oo. 

The following lemma asserts that the constant 7 defined in (fT8l) assumes different values according to 
whether G(H\, 1, W(d,t)) is in the subcritical phrase or the supercritical phase. 

Lemma 11: Let 7 be defined as (fT~8T ) . (i) If G(7i.\, 1, ^(c?, i)) is in the subcritical phase, i.e., A < 
X c (rji(d)), then 7 < 00, and 7 > with probability 1. (ii) If G(H\,l,W(d,t)) is in the supercritical 
phase, i.e., A > X c (r)i(d)), then 7 = with probability 1. 

Proof: To show (i), note that 7 < 00 follows directly from 

= inf < 1 < Mi (2g) 

m>l m 

where the last inequality is shown above in the proof for Lemma [6] 

To see why 7 is positive with probability 1, suppose the node at X disseminates a message at time 
t — 1 and consider G(H\, 1, W(d,t ))- Choose K large enough such that c\e~ C2K < |, where ci and c 2 
are the constants given in Proposition [3l Let q = I jTR+T) J ■ When m > 2(K + 1), g > 1. 

Let ^ = {(x,y) E R 2 : K + (/i - + 1) < x - z(Xq) < + 1)} for /i = 1,2,..., where ar(u) 
is the x-coordinate of node f . Since X and X m are both in C(G(TC\, 1)), there exists at least one path 
from X to X m . Moreover, since each strip S h has width 1, at least one node of C(G(Tt\, 1)) lies inside 
each Sh,. 

Let {X} 1} ,/ = 1,2,...} be the nodes of C{G{H X , 1)) which lie inside Si. Since G(H X , 1, W(d, t )) 
is in the subcritical phase, by Proposition [3l the probability that there exists a path consisting of only 
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active links from X to any Xj , I — 1, 2, is less than or equal to Cie~ C2ft: < |. In other words, with 
probability strictly greater than |, there exists at least one inactive link at time t = t on any path from 

X to x[ 1} , / = 1, 2, .... Let T« = mf z {T(X , x[ 1} )}. Let iy (<iit) 4 inf 0<d <i { Vo (d)E[Y k (d)]} > 0, then 
£[r«] > \T w{d>t) > 0. 

Let {Xj, 1 ^, Z' = 1,2,...} be the nodes of C(G(H\, 1)) which lie inside S^+i, for h > 1. By the same 
argument as above, the probability that there exists a path consisting of only active links from any node 
in Sh to any node in Sh+i is less than or equal to c\t~ c ' lK < |. In other words, with probability strictly 
greater than |, there exists at least one inactive link on any path from any node in Sh to any node in Sh+i- 
Let T {h+l ^ = infi ) ^{T(x{ ft) ) Xf +1) )}, Then Sp 1 ^ 1 )] > \T W (d,t) > 0. The path segments are illustrated 
in Figure [6j 

Since ||X — X m || > m — r — r m , when y > r + r m , any path from X to X m has at least Li(FTi)J = 9 
segments and the delay on each segment is strictly greater than |lV(<i,t) > 0. Hence, E[T 0jm ] > \qTw{d,t) 
when y > r + r m . Since both r and r m are finite with probability 1, y > r + r m holds with probability 
1 as m — >• oo. 

Since if is finite and Tw(d,t) is positive and independent of m, we have 

E[T 0rn ] 
7 = hm 

m^oo m 

> lim — -T W (d,t) 

> Hm (2(^fl) " ^ ^^(^ 



m^oo 



> (29) 

with probability 1, where we used the fact that q > 2 (k+i) ~ ^ - 

For (ii), suppose G(H\, l,W(d,t)) is in the supercritical phase. To simplify notation, let C(t) be the 
infinite component of G(H.\, 1, W(d,t)). Let t' be the first time when some node in C(t') receives X 's 
message, and let 

w\ = argminc/(Xj, Xo), and W2 = axgmind(Xj, X m ). 

ieC(t') ieC(i') 

That is, Wi and u; 2 are the nodes in the infinite component of G(TC\, 1, W(d, t')) with the smallest 
Euclidean distances to nodes X and X m , respectively. If node X is in C(t ), then t' = t and w\ = X . 
If at time t', node v is in C(t'), then u> 2 = X m . 

Since both wi and iy 2 belong to C(t'), T(wi, W2) = 0. The distances d(Xo,X u , 1 ) and c?(X W2) X m ) are 
finite with probability 1 by Lemma [16] in Appendix C. Clearly, rf(X ,X Wl ) is independent of m. By 
stationarity, d(X W2 ,X m ) is also independent of m. Hence, by the proof of Lemma [6l E[T(X. ,'K Wl )] < 
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oo, E[T(X W2 ,X m )) < oo with probability 1 for any m, and E[T(X ,X Wl )] and E[T(X W2 ,X m )) are 
independent of m. Moreover, 



< T 0tin ^ T(X , X Wl ) + T(wi, w 2 ) + T(X W2 , x r 
~~ m ~ m 

T(Xq, X Wl ) + T(X W2 , X m ) 



m 



(30) 



Ell} 1 

Hence 7 = linv^oo 1 0,mJ = with probability 1 . □ 



We are now ready to prove Theorem [5J 

Proof of Theorem \5\- Assume node u disseminates a message at time t = to- Take X u as the origin, 
and the line X U X„ as the a;-axis. By definition u,v E C(G(H\, 1)). Since node u is the origin, X u = X . 
Let m be the closest integer to x(v) — the x-axis coordinate of node X„. Now T 0jm = T(X u ,X m ). If 
X„ = X m , T{u,v) = T , m . 

Note that m — 1 < d(u, v) < m + 1, Thus, for any m > 1, we have 

,m /o 1 \ 



m + 1 d(tt, f ) m — 1 
On the other hand, if X„ 7^ X m , then X m must be adjacent to X„. This is because 1 1 (m, 0) — X„| | < | (m 
is the closest integer to x(v)) and ||(m, 0) — X m || < | (X m is the closest node to (m, 0)). Consequently, 
T , m - T(X m , X,) < T(n, u) < T , m + T(X m , X,). Thus, for any m > 1, we have 

7p,m - r(x m , x^) ^ t(u,v) ^ r 0[m + T(x m , x B ) ^ 32 ^ 

m + 1 <i(w, f ) m — 1 

Since X m is adjacent to X , T(X m , X„) < 00 with probability 1. Therefore, in both cases, by Lemma[6] 
and a typical e-S argument (see Appendix D), we have for any e > 0, 5 > 0, there exists d < 00, such 
that if d(u, v) > do, then 



Pr 



T(u,v) 

7 



< e > 1 - 5. (33) 



d(u, v) 

When G(H\, 1) is in the subcritical phase, by Lemma ITTl we have < 7 < 00 with probability 1. 

On the other hand, when G(H\, 1) is in the supercritical phase, by Lemma ITTl we have 7 = with 
probability 1. Then, by a typical e-S argument (see Appendix E), we have for any e > 0, 5 > 0, there 
exists d < 00, such that if d{u, v) > d then 

Pr ri^4 <£ ) >1 -i. 

d(u,v) ) 

□ 
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C. Effects of Propagation Delay 

Up to this point, we have ignored propagation delays. We now take this type of delay into account. 
Suppose the propagation delay is < r < oo for any link, independent of the link length. We assume 
the following mechanism is used for a transmission from node i to node j: (i) a packet is successfully 
received by node j if the length of the active period on link during which the packet is being 

transmitted, is greater than or equal to r; (ii) node i retransmits a packet to node j until the packet is 
successfully received by j. 

Note that due to the Markovian nature of the link state processes {Wij(dij,t)}, at the instant when a 
packet arrives at node i, the residual active time for link has the same distribution as Z(dij). Thus 
without loss of generality, we assume that node i initiates transmission on link at time 0. If link 
is on at time with Zi(d) > r, then the transmission delay T[Ad) on is r. However, if link 
is on at time with Z\(d) < r, or if is off at time t = 0, then the delay on is less 
straightforward to calculate. In this case, we need to capture the behavior of retransmissions. Let 

K{d) = argmm{Z k (d) > r}. (34) 
fc>i 

Then, K(d) is a stopping time for the sequence {Z k (d), k > 1}. Now we have 

K-l 

T[ j = J2( Y * + Z *) + Y K + r > W(d,0) = 0, 

k\ (35) 

^ = E^ + Z *) + r ' W(d,0) = l, 

i=i 

where we abbreviate T^{d), K(d), Yi(d) and Zi(d) as T[j, K, Yi and Z, t , respectively. 
Let 

T T {u,v) = 7-(X u , X„) 4 inf \ W . 06) 

where T[j(dij) is given by (1351) . Then, T T (u,v) is the message delay on the path from u to v with the 
smallest delay, including propagation delays. 

Corollary 12: Given G(H\, 1, W(d, t)) with A > A c and propagation delay < r < oo, there exists a 
constant 7(7-) < 00 with j(t) > r (with probability 1), such that for any u,v E C(G(H\, 1)), and any 
e > 0, 5 > 0, there exists d < 00 such that for any u, v with d(u, v) > d , 

T T {u,v) 



Pr 



d{u,v) 7(r) 



< e > 1 - 5. (37) 
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Moreover, when G(H\,l,W(d,t)) is in the subcritical phase, as r — > 0, j(t) — > 7 with probability 
1, where 7 is defined in Theorem \5\ When G(H\,l 1 W(d,t)) is in the supercritical phase, as r — > 0, 
7(r) ^ with probability 1. 

To prove this corollary, we need the following two lemmas. 

Lemma 13: Given any < r < 00, for all < d < 1, the expected delay on each link is positive 
and finite, i.e., 

< E[T^\ < 00. (38) 

Proof: By (1331) , we have 

E[I%] = E[E[T%\W{d,ti)]\ 

= rioEimW&O) = 0] +r ]l E[T[ j \W(d,0) = 1] 

= VoE 



K-l 



+V1E 



J^iXi + Zi) +Yk + r\Zi < r,i = 1, X - 1 
=1 

A'-l 

^(Ki + Z t ) + r|Z, < r, 2 = 1, K - 1 



i=l 



T + r] E[Y K ]+E 



K-l 



J2{Yi + Zi)\Zi<T,i = l,...,K-l 



i=l 



< E[K)r + Vo E[Y K ] + (E[K] - l)^], 



(39) 



where in the last equality, we used the fact that Yi and Zj are i.i.d. and Zj < r for z = 1, 2, ...if — 1, as 
well as Wald's Equality for stopping time K. 

Since < r < 00, < T] < 1, and < < 00, in order to show < E[T^] < 00, it suffices to 

show 1 < E[K] < 00. By definition, K > 1 so that i£[iT] > 1. Thus, we need only to show E[K] < 00. 
For any k > 1, Pr(fsT = jfe) = Pr(Z a < r, ...,Z fc _i < r, Z fe > r) = F z (r)* _1 (l - ^(r)), where 
F z (-) = Pr(Z < r). Then 



00 

= kF z {r) k -\l - F z (r)) = — 

k=l 



Fz(r) 



Therefore, we have E[K] < 00. 



(40) 
□ 



Lemma 14: Given G(TC\, 1, W(d, t)) with A > A c and no propagation delay, let L , m be the path from 
Xo to X m that attains T 0jm and has the smallest number of links (in case there exist multiple paths 
attaining T 0jm ). Then |L 0iTO | < 00 with probability 1 for each m, where |L 0jm | is the number of links 
along L 0>m . 
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Proof: By the proof of Lemma [91 we have E[T 0jTn ] < oo. We can express E[T 0jTn ] as 

E[T 0:Tn ] = E[E[T 0)m \\Lo >m \)), 

where 

E[T 0;m \\L 0;m \] = ^\d)E[Y^\d)} > \Lo >m \V wm , 
i=i 

where i]q \cI) and E[Y^ l \d)] are the stationary probability of the inactive state, and the expected inactive 
period of the i-th link with length d on L 0m respectively, and T w ^ t ) = in£ 0<d < 1 {rio(d)E[Yf,(d)}} > 0. 
Thus, we have 

E[\L 0tm \]T W ( d ^ < oo. 

This implies i?[|L ,m|] < oo, which further implies |Lo,m| < 00 w i m probability 1. □ 

Proof of Corollary \U\ Let T£ m = T T (X,,X m ), for ||X, - X m || < 00, < I < m, where % is denned 
as in (f!6l) . 

Clearly, the relationship TJ" m < TJ" z + T[ m still holds for any < I < m. Hence, condition (i) of 
Theorem |7J holds. Since the propagation delay does not affect the stationarity of the geometric structure 
of the network, conditions (ii) and (iii) of Theorem |7J also hold. 

By the same argument as that in the proof of Lemma[9l we have < 00, where \L\ = |L(X , X m )| 

and L(X , X m ) is the shortest path from X to X m . Let Tq be the delay on this path. Then, 

\L\ 

E[T^\\L\] =Y J E[Tl{di)\ < WA^t,, 
i=i 

where T[(di) is the delay on the i-th link with length di on the path L(X , X m ), as given by (|35l) . and 

Aw^(d,t) — sup 0<d<1 E[T[(di)} < 00. By Lemma [HI we have < E[T[(di)} < 00 for all < di < 1, so 
that hw T {d,t) < 00. Hence 

E[T^} = E[E[T^\\L\}} < E[\L\]A W r m < 00, 

which implies E[TQ m ] < 00. This ensures that condition (iv) of Theorem |7j holds. 

Furthermore, the propagation delay does not affect the strong mixing property of {TJ r m , < I < m}. 

a E[T T ] 

Therefore the result of Lemma [6] holds for {Tf m , < I < m}. Let 7(7-) = lim™-^ — then 7(1-) = 
inf m >! ^=1, and 

rpj- 

lim — — = j(t) with probability 1. (41) 

TO— >oo m 
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Then applying the same proof for Theorem [51 we can show that for any e > 0, 5 > 0, there exists 
do < oo, such that if d(u, v) > d , then 

T T (u,v) 



Pr 



7(t) 



< e > 1 - 5. 



d(u, v) 

To see why j(t) < oo and 7(r) > r with probability 1, first note that 

7 t = ™f < ^il < oo. (42) 

m>l m 

Moreover, since the shortest path between nodes X an d X m has at least |_| |X — X m | |J > [m — r Q — r m \ 
links, T T m > t [m — r — r m \ . Since r and r m are both finite with probability 1 and independent of m, 
we have 7(7-) > r with probability 1. 

In the following, we show that as r — > 0, 7(r) — ► 7 with probability 1 when G(H\, 1) is in the 
subcritical phase, and 7(7-) — > with probability 1 when G(H.\, 1) is in the supercritical phase. Observe 
that 

1-^0, m| 

T , m < < ]T 27(d ( ), 

where L 0>m is defined in Lemma [I4l and T[(di) is the delay on the i-th link with length dj along L , m , 
as given by (|35l) . From Lemma O we have |L 0>m | < 00 with probability 1. Thus with probability I, 

\Lo,m\ 

E[T , m ] < E[TS,J < E[T[(di)]. 
By (El and E[T , m ] = E&i™ 1 77o(^)£[W)] we have 

1-^0, m| 

£[T , m ] < E[rj; m ] < £[T , m ] + ILo^I^^t + ( E l R ] - (43) 

i=l 

with probability 1. From (l40l) . we know that as r — > 0, E'fif] — > 1. Therefore, as r — > 0, we have 
|L ,m|£[#]T + El=i m '(£[^] - l)^[^fc(^)] -> with probability 1. This, combined with @3 implies 
lim r _^o E[Tq m ] = E[T 0im ] with probability 1. Therefore, 

lim7(r) = lim lim — - 

hm lim 



m^oo r^0 777, 

= hm 

m^oo 777, 

= 7, (44) 

with probability 1, where the interchanging of limitation operations is justified by E[TQ m ] < oo. Conse- 
quently, as r — ► 0, 7(r) — ► 7 with probability 1 when G(TC\, 1) is in the subcritical phase. Since 7^0 
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Fig. 7. Delay performance of information dissemination in wireless networks with dynamic unreliable links (A = 1.75): (a) E[Ti(d)] = 0.5 
and E[T (d)} = 2 for any < d < 1; (b) E[Ti(d)] = 2.5 and E[T (d)] = 0.5 for any < d < 1. 

with probability 1 if G(H\, 1) is in the supercritical phase, we have 7(7-) — > with probability 1 in this 
case. □ 



An interesting observation of this corollary is when the propagation delay is large, the message delay 
cannot be improved too much by transforming the network from the subcritical phase to the supercritical 
phase. However, as the propagation delay becomes negligible, the message delay scales almost sub-linearly 
(7(1") ~ 0) when the network is in the supercritical phase, while the delay scales linearly (7(7") ~ 7) 
when the network is in the subcritical phase. 

V. Numerical Experiments 

In this section, we present some simulation results. Figure UM\ show simulation results of the information 
dissemination delay performance in large-scale wireless networks with dynamic unreliable links. 

In Figure U\ the lengths of the active and inactive periods have exponential distributions independent 
of d — the length of the link. In Figure [H the lengths of the active and inactive periods have exponential 
distributions depending on d. In all of these scenarios, it can be seen that when the resulting dynamic 
network is in the subcritical phase, converges to a non-zero value as d(u,v) — > 00. The limit 

depends on the density of G(H\, 1) and the distributions and expected values of the active and inactive 



periods. When the resulting dynamic network is in the supercritical phase, 



T(u,v) 
d(u,v) 



converges to zero as 



d(u, v) — > 00. 

To see how propagation delays affect the message delay, and to verify the results of Corollary [121 we 
illustrate simulation results in Figure [9l where 71(d) and T (d) have exponential distributions independent 
of d. 
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Fig. 8. Delay performance of information dissemination in wireless networks with dynamic unreliable links (A = 1.875): (a) E[Ti(d)] = 0.5 
and E[T (d)} = 1.5d + 1 for any < d < 1; (b) £pi(d)] = 2 and E[T (d)} = 0.5d + 0.5 for any < d < 1. 
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Fig. 9. Delay performance of information dissemination in wireless networks with dynamic unreliable links (A = 1.875) and propagation 
delay r = 1: (a) E[Ti(d)] = 1 and E[T (d)} = 8 for any < d < 1; (b) E[Ti(d)} = 1 and E[T (d)} = 2 for any < d < 1. 



VI. Conclusions 

In this paper, we studied percolation-based connectivity and information dissemination latency in large- 
scale wireless networks with unreliable links. We first studied static models, where each link of the network 
is functional (or active) with some probability, independently of all other links. We then studied wireless 
networks with dynamic unreliable links, where each link is active or inactive according to Markov on-off 
processes. We showed that a phase transition exists in such dynamic networks, and the critical density 
for this model is the same as the corresponding one for static networks (under some mild conditions). We 
further investigated the delay performance in such networks by modelling the problem as a first passage 
percolation process on random geometric graphs. We showed that without propagation delay, the delay of 
information dissemination scales linearly with the Euclidean distance between the sender and the receiver 
when the resulting network is in the subcritical phase, and the delay scales sub-linearly with the distance 
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if the resulting network is in the supercritical phase. We further showed that when propagation delay 
is taken into account, the delay of information dissemination always scales linearly with the Euclidean 
distance between the sender and the receiver. 

Appendix A 

Proof of Proposition]^ Let B be a bounded box containing the origin, and let W(B) be the union of 
components that have some node(s) of G(H\, 1, p e (-)) inside box B. Precisely, W(B) = {component W G 
G(H x ,l, Pe (-)):3we W'^eB}. 

Consider the following two events: 

E = {d(W(B)) > h}, and F = {all nodes of G(H X , l,p e (-)) in side B belong to W }. 

Clearly, events E and F are both increasing events. By the FKG inequality, we have Pi(E n F) > 
Pr(E) Pr(F). Thus, 

Pr(d(W ) > h) > Pr(E n F) 
> Pr(E)Pr(F) 

= Pr(F) Pr (d(W(B)) > h), (45) 

where Pr(F) > since B is bounded. By (|45l) . we have 

E[d(W )] 



E[d(W(B)\ < 



Pr(F) 

Therefore, when A < A c (p e (-)), we have E[d(W )] < oo and thus E[d(W(B)] < oo. 

To prove the Proposition, it is sufficient to show Pr(£> <~~> B{h) c ) < cie~ C2h , where {B B{h) c } 
denotes the event that some node(s) inside B and some nodes in B{h) c are connected. 

We partition the space as the union of B(i,j) = (i — |, % + |] x (j — \, j + |] , where (i, j) E Z 2 . Since 
E[d(W(B(0, 0))] < oo, d(W(B (0,0)) < oo with probability 1. Then we can choose M sufficiently large 
so that E[H M ] < |, where H M is the number of boxes B(i,j) outside B(M) = [— M, M} 2 intersecting 
W(B(0,0)). 

Now choose L large enough so that the set U m (i j)>l-i B(i,j) is disjoint from B(M), where m(i,j) = 
max.{\i\,\j\}. Choose h sufficient large so that U m uj) <L B(i,j) C B(h). Observe that if {5(0,0) <-~» 
B(h) c } occurs, then there exists (i,j) with m(i,j) = L for which {5(0,0) <^ D(i,j)} and {B(i,j) 
B(h) c } occur disjointlyj^, where D(i,j) — ')=l— i m{i—i' j— j')=i Bty'ij')- This is illustrated in Fig- 

6 Let U be a bounded Borel set in R 2 . For any realization G G G(H X , l,p e (-))> let G u = (K,,-E u ), where V u = {v : v € G P U} 
and E u = {(u, u) : it, w £ K}- Define [G u ] = {G" £ G(i?x, l,Pe(-)) : 3G" C G' s.t. G^ = G u }. We say that an increasing event A 
is an event on U if Ia(G) = 1 and G' € [G u ] imply that Ja(G') = 1. A rational rectangle is an open 2-dimensional box with rational 
coordinates. Let A and B be two increasing events on U, and W\ and W2 be two disjoint sets that are finite unions of rational rectangles. 
For G £ G(H\, l,p e (-)), if ^(G^) = 1 where G' Wl G [G Wl ], and 7 s (Gw 2 ) = 1 where G' W2 G [Gw 2 ], then we say that A and B 
occur disjointly. We use AOB to denote the event that A and _B occur disjointly. For details, please refer to [3], [4]. 
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Fig. 10. B(h), B(M), L, B(i,j) and D(i,j). 

ureM 

Let {£(0,0) «~* D(i,j)DB(i,j) ^ B (hf) denote the event that {5(0,0) D(i, j)} and {B(i,j) 
B{h) c } occur disjointly. It then follows from the BK inequality [3], [4] that 

Pr(£(0,0) B(h) c ) < Pr(5(0,0) > D(i, j)DB(i, j) B(h) c ) 

(i,j):m(i,j)=L 

< max Pr(S(i, j) B{h) c ) V Pr(S(0, 0) +~* D(i, j)) 

(i,j):m(i,j)=L ' * 

(i,j):m(i,j)=L 

= r . max Pi (B(i,j)<^B (h) c ) V E[I {Bm ^ D{iJ)} } 

(i,j):m{i,j)=L f 

(i,j):m(i,j)=L 



max Pr(B(i,j) B(h) c )E 

(i,j):m(i,j)=L 



! {B(0,0)^D(i,j)} 



(i,j):m(i,j)=L 



< max Pr(5(i,j) £(fc) c )3£[# M ] 

(i,j):m(i,j)=L 



(46) 



where the last inequality follows from the fact the each box B(i',f) can be contained in at most 3 
D(i,j)'s. 

It follows that 



Pr(£(0,0) *~* B(h) c ) < - max Pr(B(i,j) B(h) c 

2 (i,j):m(i,j)=L 



(47) 



To bound the right hand side of (07]), choose a sufficiently large h such that Um(i'-ij'-i)=L,m(i,i)= 
£?(/i). The same argument as above shows that for all with m(i,j) = L, 



Pv(B(i,j)^B(h) c )<- max Pr(fl(i', /) ^ B(h) c ). 

Z (i',j'):m{i'-i,3'-j)=L 



(48) 
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Fig. 11. An example of a circuit surrounding the origin in lattice £' 

Repeating this argument leads to the desired conclusion. □ 



Appendix B 

The following lemma is similar to the one used in [3], [9], [11]. For completeness, we provide the 
proof here. 

Lemma 15: Given a square lattice £', suppose that the origin is located at the center of one square. 
Let the number of circuits^ surrounding the origin with length 2m be 7 (2m), where m > 2 is an integer, 
then we have 

T(2m) < ^(m-l)3 2m (49) 

Proof: In Figure [TT] an example of a circuit that surrounds the origin is illustrated. First note that the 
length of such a circuit must be even. This is because there is a one-to-one correspondence between each 
pair of edges above and below the line y = 0, and similarly for each pair of edges at the left and right of the 
line x = 0. Furthermore, the rightmost edge can be chosen only from the lines li : x — i— |, i — 1, m— 1. 
Hence the number of possibilities for this edge is at most m — 1. Because this edge is the rightmost edge, 
each of the two edges adjacent to it has two choices for its direction. For all the other edges, each one 
has at most three choices for its direction. Therefore the number of total choices for all the other edges 
is at most 3 2m ~ 3 . Consequently, the number of circuits that surround the origin and have length 2m must 
be less or equal to (m — l)2 2 3 2m ~ 3 , and hence we have (l49l) . □ 



7 A circuit in a lattice £' is a closed path with no repeated vertices in 
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Appendix C 

Lemma 16: Suppose G(H\, l,p e (-)) is in the supercritical phase, i.e, A > A c (p e (-)). Let v C(G(H\, l,p e 
and define 

w = argmin d(i,v), 

ieC{G(Hx,l, P e(-))) 

i.e., w is the node in the infinite component of G(H\, l,p e (-)) with the smallest Euclidean distances to 
node v. Then, d(w,v) < oo with probability 1. 

The idea behind the proof for this lemma is similar to that for the proof for Lemma [8] The difference 
is that the probability of a good event is now defined with respect to G(H\, l,p e (-)) instead of G(H\, 1). 

Given G(H\, l,p e (-)) with A > A c (p e (-)), as in the proof for Lemma [H we consider a mapping between 
G(H\, l,p e (-)) and a square lattice £ = d- Z 2 , where d is the edge length. The vertices of C are located 
at (d x i, d x j) where G Z 2 . For each horizontal edge a, let the two end vertices be (d x a x ,dx a y ) 
and {d x a x + d, d x a^). 

As in the proof for Lemma [H define event A a (d,p e (-)) for edge a in £ as the set of outcomes for which 
the following condition holds: The rectangle R a = [a x d— |, a x d+^} x [a y d—^, a y d+f\ is crossed from left 
to right by a connected component in G(H\, 1, p e (-)). Define event A! a {d, p e (-)) for edge a in £ as the set of 
outcomes for which the following condition holds: The rectangle R a = [a x d—^, a x d+^] x [a y d—^, a y d+^\ 
is crossed from top to bottom by a connected component in G(H.\, l,p e (-)). 

Let 

p g (d,p e {-))±Pr{A a {d, Pe (-))l and p' g (d,p e {-)) ± Pr{A' a {d,p e (-))). (50) 

Define A a (d,p e (-)) and A' a (d,p e (-)) similarly for all vertical edges by rotating the rectangle by 90°. 

Define a vacant component V in IR 2 with respect to (w.r.t.) G(H\, l,p e (-)) to be a region V CM 2 such 
that V H G(H\, l,p e ( - )) = ^ (i- e -' no no( i e or an Y P ar t of a link of G(H\, l,p e (-)) is contained in V), and 
such that there exists no other region U C M 2 satisfying Vc C/ and £7 fl G(TC\, l,p e (-)) = 0. 

Definition 3: For G(H\, l,p e (-)), let V be the vacant component in R 2 w.r.t. G(H\, l,p e (-)) containing 
0. Let 

A:(p e (0) = sup{A : Pr(rf(V ) = oo) > 0}. (51) 

Similarly we can define the vacant component Vq containing the origin in R 2 w.r.t. G(H\, 1), and 
A* = sup{A : Pr(d(V ') = oo) > 0}. It is known that A* = A c (Chapter 4 in [4]). Since G(H\, l,p e (-)) is 
a subgraph of G(H\, 1), it is clear that A*(p e (-)) > A*. 
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Proposition 17: Let ip*(p e (-)) = Pr(3 vacant component Kcl 2 w.r.t. G(H\,l,p e (-)) ■ d(V) = oo). 
Then 

**wo) ={J; <w 

Proof: First assume A < A*(p e (-)). The graph G(H\, l,p e (-)) is obtained by placing a link between 
two nodes % and j with probability p e (') when ||xj — Xj|| < 1. The event {3 vacant component V C 
IR 2 w.r.t. G(TC\,l,p e (-)) : d(V) = oo} does not depend on the existence of any finite collection of 
those links. By Kolmogorov's zero-one law [3], [22], ip*(p e (-)) assumes the values and 1 only. Since 

Pr(d(Vb) = oo) > 0, then 

riPe(-)) > Pr(rf(Vo) = oo) > 0, 

so that ip*(p e (-)) = 1 by Kolmogorov's zero-one law. 

On the other hand, if A > A*(p e (-)) > A c , with probability 1, there is no vacant component with infinite 
diameter in IR 2 w.r.t. G(H X , 1) (Chapter 4 in [4]). Since Pr(d(V r ) = oo) = 0, we have 

r(Pe(-))< £ Pr(d(Vy = °o) = 0, 

xeQ 2 

where we used the fact that Q 2 is dense and any infinite vacant component is open so that any infinite 
component contains at least one x e Q 2 . □ 

Given the mapping between G(H\, l,p e (-)) and C, define event A*(d,p e (-)) for edge a in £ as the set of 
outcomes for which the following condition holds: the rectangle R a = [a x d— |, a x d+^f] x [a y d—^, a y d+^\ 
is crossed from left to right by a vacant component in IR 2 w.r.t. G(H\, l,p e (-)). Define event A^(d,p e (-)) 
for edge a in £ as the set of outcomes for which the following condition holds: the rectangle R a = 
[a x d — |, a x d + ^] x [a y d — |, a y d + |] is crossed from top to bottom by a vacant component in R 2 w.r.t. 
G(H x ,l, Pe (-)). 

Let 

p;(d :Pe (-))±Pr(A:(d,p e (-))), and p^{d,p e {-)) 4 Pr(^(d lPe (.))). (53) 

Define A*(c?,p e (-)) and A*'(d,p e (-)) similarly for all vertical edges by rotating the rectangle by 90°. 
Figure [121 illustrates A*'(cf,p e (-)). 

We now define another critical density with respect to G(H\, l,p e (-)). 

Definition 4: Given G(H\, l,p e (-)), let 

A|(p e (-)) = sup{A : lim sup p*g(d, Pe {-)) > 0}. (54) 
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Proposition 18: For G(H\, l,p e (-)), we have 



x c ( Pe (-)) = K( Pe (-)) = y s (Pe(-))- (55) 

Proof: To show (O, it is sufficient to show (i) A c (p e (-)) < A*(p e (-)), (ii) A*(p e (-)) < A£(p e (-)), and 

(iii) AJ(p (.)) < A c (p e (0). 

To show (i) A c (p e ( - )) ^ A*(p e ( - ))' l et A < A c (p e (0)- Then G(H\,l,p e (-)) is in the subcritical phase. 
Let Bi(i) = (0, 2i) + B(l) where 5(1) = [-1, l] 2 for z = 0, 1, 2.... Observe that the existence of a left to 
right crossing in rectangle [0,3 fc ] x [0, 3 fc+1 ] by a component W of G(H\, l,p e (-)) implies the existence 
of a component VT" of G(H\, l,p e (-)) starting from |Ji=o 2 ^iv) (i- e -> me fi rst n °d e m W m the x-axis 

r 3*+ 1 i 

direction is inside |J i=0 2 Bi(i)) with diameter greater than or equal to 3 fc — 2. Hence, we have for any 

k > 1, 



p'(d = 2-3 fc ,p e (-)) < Pr 



Q K^(5 1 (.)))>3 fc -2} 



i=0 



/ 



gk + 1. 

2 



< [J Pr(d(W(Bi(i))) > 3 fc - 2) 

i=0 

•gjb+n 



< 



< 



+ 1 Pr(d(W(B(l))) > T -2) 



fc-i 



2 J Pr(d(W(5(l))) > 3 fc -2) 
2 J Pr(d(W(B(l))) > 3^), 



(56) 



where W(B 1 (i)) is the union of components of G(H\, l,p e (-)) that have some node(s) inside box B 1 (i). 
Precisely, W(B x (i)) = {component W of G(H X , l,p e (-)) : 3w e W',*^, e 

Since A < A c (p e (-)) = \£>(p e (-)), E[d(Wo)) < oo. By the same argument used in the proof for 
Proposition we have E[d(W(B(l))) < oo. 
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Fig. 13. A vertical crossing t k of [0,3 fc ] x [0,3 fc+1 ] and a horizontal crossing l k of [0, 3 fc+2 ] x [0, 3 fc+1 ] must intersect. 

Let P k = Pi(d(W(B(l))) > k). Then P k is non-increasing in k, and thus we have 

oo oo , q \ 

Y,p' 9 (d = 2-3 k ,p e (-)) < ^ -3^ + 2 P 3 - 

fc=l k=l ^ ' 

; — n \ / 



k=0 k=0 

^ | ( p i + 3 E 3 *" lp ^ J + 2£[d(W(B(l)))] 

< |(P! + 3E[d(W(B(l)))}) + 2E[d(W(B(l)))} 

< oo. 



(57) 



Note that p' g (d = 2 • 3 k ,p e (-)) + p* g (d = 2 ■ 3 k ,p e (-)) = 1 for all k > 1. Hence by the Borel-Cantelli 
Lemma, we have 

Pr(3 vacant top to bottom crossing t k in [0, 3 fc ] x [0, 3 fc+1 ] for all suffcient large k) = 1. 
Rotational invariance implies that 

Pr(3 vacant left to right crossing l k in [0,3 fc+2 ] x [0,3 fc+1 ] for all suffcient large k) = 1. 

As illustrated in Figure \\3\ a vertical crossing t k of [0, 3 fc ] x [0, 3 fc+1 ] and a horizontal crossing l k of 
[0, 3 k+2 ] x [0, 3 fc+1 ] must intersect. Also, t k+1 of [0, 3 fc+1 ] x [0, 3 k+2 ] and l k must intersect. Thus the union 
of vacant crossings {t k } and {l k } combines to give an infinite vacant component in the first quadrant. 
Therefore, by Proposition [13 A < A*(p e (-)), and A c (p e (-)) < A*(p e (-))- 

We now show (ii) A*(p e (-)) < AJ(p e (-)). Let A > X* s (p e (-)). Then limsup d _ >00 p*'(tZ,j) e (-)) = 0, and 
hence limsup d ^ 00 p g (d,p e (-)) = 1. Then there exists 5 > such that there are infinitely many d[, d' 2 ,... 
satisfying p g (d' i ,p e (-)) > 8 for i — 1, 2, .... Now choose d\ = d[ and = min{c^ : d'- > 3d-}. Then by 
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the same argument used in the proof for Lemma |9l we can construct infinitely many annuli around the 
origin, each annulus having edge length d\ and containing a circuit with a probability larger than 5. Then, 
by the Borel-Cantelli Lemma, with probability 1, there exist infinitely many circuits surrounding the origin 
and hence d(V ) is finite with probability 1. This implies that A > A*(p e (-)), and thus X* s (p e (-)) > X*(p e (-)). 

Finally, (iii) X* s (p e (-)) < A c (p e (-)) can be shown by the same argument as that for the proof of Theorem 
4.3 and Theorem 4.4 in [4]. □ 



Proof of Lemma I76l- If G(H\,l,p e (-)) is in the supercritical phase, A > X c (p e (-)) = A*(p e (-)) = 
X* s (p e (-)). Thus, limsup d ^ 00 p*'(ci,p e (-)) = and limsup d ^ 00 p £ ,(<i,p e (-)) = 1. Then by the same methods 
used in the proof for Lemma [81 we can show Lemma [16] □ 



Appendix D 

Since T(X m ,X t; ) < oo with probability 1, for any < 5i < 5, there exists M < oo such that 

Pr(T(X m ,X u )<M)>l-5 1 . 

Then for any e > 0, 



Pr 



T(u,v) 



d(u, v) 



7 



■> = p, r f< « <i+f , 

/ T(u v) 

> Pr 7 - e < < 7 + e|T(X m , X v ) < M ) Pr(T(X m , X„) < M) 



> 


Pr 


v 7 " 


e < 




Pr ( 


:(7- 


e + 


> 


Pr ( 


'(7 + 


e + 



T(u 


v) 


d(u, 


v) 


T(u 


V) 


d(u, 


v) 


T(u 


v) 


d(u, 


v) 


^0,m 





> Pr(7-e<3^T<7 + ^(X m ,X u )<Ml (1-6,) 

< 7 + e ] (1-Si) 



m + 1 m — 1 



Since linv^oo = 7 with probability 1, for 8 2 = 1 — there exists m < 00 such that for any 

m > m , 



2 m 2 



If 7-f <^<7 + f,then 



T , m < (7 + |) m < m( 7 + e) - (M + 7 + e), 
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and 



T , m > (7-|)m>m(7-e) + (M + 7 + e) 



Hence, for any m > max{m , 2 ( M +7+^) we j^ye 

Pr ((7 + e + M) + (7 - e)m < T , m < 771(7 + e) - (M + 7 + e)) > 1 - 5 2 . 
Moreover, since m > d(u, v)—l, if d(u, v) > d = max{m , 2 ^ Af+7+ ^ }+l, we have m > max{m , 2( - M+7+£ - > }, 



so that 



Pr 



T(w,u) 



7 



<e > (l-(Ji)(l-(J 2 ) = 1-5. 



Appendix E 

Let e>0, 0<5<lbe given. When G(H\,l,W(d,t)) is in the supercritical phase, 7 = with 
probability 1. Thus, there exists < e\ < e and < 81 < 8 such that 

Pr(7 < ei) > 1 - 81. 

Let e 2 = e — ei, and 5 2 = 1 — yz^"- F rom Appendix D, we know that for e 2 and 8 2 , there exist d < 00 
such that when <i(w, f ) > d , 

T(u, v) 



Pr 7 - e 2 < 



d(u, v) 



< 7 + e 2 > l-5 2 . 



Then for the given e, when d(u,v) > d , we have 

„ (T(u,v) \ 
Pr ( .' ' / < e > Pr 



d(u, v) 



> Pr 



T(u, v) 


d(u, 


v) 


T(u 


v) 


d(u, 


v) 


T(u 


v) 


d(u, 


v) 



< e|7 + e 2 <ej Pr(7 + e 2 < e) 

< e|7 + e 2 < e J (1 -81) 



< 7 + e 2 (1 - 81) 



> Pr 

> (1-8 2 )(1-8 1 ) 
= 1-8. 
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